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In order to improve simulations with domain wall fermions (DWFs), it has been suggested to project out a 
number of low-lying eigenvalues of the 4-dimensional Dirac operator that generates the transfer matrix of DWF. 
We investigate how this projection method affects chiral properties of quenched DWF. In particular, we study 
the behaviour of the residual mass as a function of the size of the extra dimension. 



1. Introduction 

Domain wall fermions (DWFs) preserve chi- 
ral symmetry [|l|,||J|] when the lattice size in the 
5th direction, TV^, is taken to infinity. A mea- 
sure of chiral symmetry breaking is the residual 
mass mros- Even though the restoration of chi- 
ral symmetry is expected to be exponentially fast 
in in practice TOres decreases very slowly as 
first shown by CP-PACS [|,|. This is due to the 
existence of very small eigenvalues of the trans- 
fer matrix along the 5th direction. At large Ns 
these eigenvalues determine the rate of the expo- 
nential decay. It is clear that any improvement 
of the chiral properties of DWF has to come from 
eliminating these low-lying eigenvalues. 

One idea is to improve the gauge actions 
11111,^,01. However, besides the potential difficul- 
ties with unitarity violations and the sampling 
of topological charge sectors, this method docs 
not solve completely the problem. In particular 
with the Iwasaki gauge action, the convergence 
rate also becomes slow at large TV^ (|] (see Fig. 
H): very small eigenvalues of the transfer matrix, 
even though less frequently, also appear in this 
case. Using the DBW2 action seems to be much 
better in this respect Q , but it is unclear whether 
these small eigenvalues could eventually appear 
there, too, leading to similar problems. 

Another method to eliminate the small eigen- 
values is to project them out of the transfer ma- 
trix H,^. In this contribution, we investigate the 
projection method based on Ref.|p|], where the 
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projection is performed on the transfer matrix it- 
self. (In Ref.||], the projection is operated only 
on a boundary term.) The aim of this paper is to 
investigate the effects of the projection method on 
the residual mass in the quenched configurations. 

2. Domain wall fermions 

Details of our notation can be found in Ref. , 
and we give only a brief summary below. 
The 5D domain wall operator is defined as 



2? = ^ {75 (9: 



(1) 



where the operator M. is obtained from the stan- 
dard 4D Wilson-Dirac operator with negative 
mass (DW mass). Choosing open boundary con- 
ditions in the 5th dimension, chiral modes with 
opposite chiralitics arc localized on the 4D bound- 
ary plane at s = 1 and N^. The action of "quark 
fields" , which are constructed from the boundary 
fermions, is related to an effective 4D operator 
DjVsj which satisfies 



A 

aD = lim aZ? w =1 , 



where 



A 



(2) 



(3) 



2 + asM 

The D of cq. (||) satisfies the Ginsparg-Wilson re- 
lation and thus an exact lattice chiral symmetry. 

In realistic simulations, for finite Ns, chiral 
symmetry is explicitly broken to a certain extent 
that can be quantified by the values of the so- 
called residual mass, which measures the anoma- 
lous term in the axial Ward-Takahashi identity 
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0: 2J^q{x) EE V^A^(x) - 2mfP{x), where 
is the axial current and P is the pseudoscalar 
density. The size of this extra chiral symmetry- 
breaking term can be described by the residual 
mass 



_ Ex(^5,ft,x)Pft,x)) 



(4) 



2.1. Improvement of domain wall fermion 

We employ a method to project out the small 
eigenvalue of A''' A The improved operator Al, 
instead of A4, satisfies the following relations in 
order that eqs.(||) and (||) hold; 



det(2 + a5Al) ^ 



(5) 



This is used to construct A4 in such a way that 
the very low eigenvalues of A disappear, while 
keeping D invariant. Following Ref. the new 
operator A4 satisfying eq.(|[) is given by the fol- 
lowing expression; 

A = '"^"^ = A + - afe)75Wfe ® vl , (6) 

2 + a^M ^ 

where Vk is the eigenvector of 75A, 

75Aufc = akVk , Js^vk = cikVk , {vk,vi) = Sm , (7) 

and r is the number of eigenvalues projected out. 
Therefore an improved DWF operator can be ob- 
tained after substituting A4 in eq. (§) with M 
given by 

r 

a^M = a^M - ^ XkiWk ® w/75 , (8) 

where {X^^)ki = 24((Sfe - afc)"^ + (vfc, w;) and 
Wfe = (2 + a5M)-f5Vk- 

3. Eigenvalues of matrix A^ A 

The eigenvalues of A^ A can be calculated by 
a straightforward generalization |^ of the Ritz 
functional method in Ref . [p^ . We have cal- 
culated the convergence rate from eigenvalues 
of A'^ A 1^,0 and the distributions of low-lying 
eigenvalues for various gauge actions pi|. 



4. Simulation of improved DWF 

The values ak can be chosen freely as long as 
Sfe > ak- This can be used to reduce the value of 
the residual mass. Here we use ctk — 2sgn(Q;A;)|a/| 
with (max{fc}; /) = (3; 10) and (10; 10) in the pla- 
quctte gauge action. The notation in the figures 
is Aproj 2A/ = 2|a/|, then = sgn(Q;fc)Aproj- 

For the Wilson gauge action we have chosen the 
gauge couphng to be /? = 6.0 (a~^ ~ 2 GeV) and 
the DW mass to be toq = 1.8. The lattice size is 
12-^ X 24 X Ns- This size is a little smaller than 
the ones used in simulations of the CP-PACS and 
RBC collaborations {16^ x 32 x Ns); the value of 
the lattice spacing is, however, the same. The 
bare quark mass is always fixed to my = 0.02. 

In Fig. I], the quantity ^'^pp-^^ without and with 
projection for Ns = 32 and 48 is plotted. The 
number of projected modes is 0, 3 and 10. Note 
the large fluctuations of the residual mass when 
no projection is performed. As the number of 
eigenvalues projected out is increased, the deter- 
mination of the residual mass is much more stable 
and cleaner. 



<J5q.P>/<PP> vs. time 

=6.0, 12\24xN„ m0=1.8, mf=0.02 

C Ns=32 (No ProjO 

Ns=32 (Proj.=3, Xp,„|=2i,„) 

►■Ns=32(Proi.=10, Vi=2^,.) 
*Ns=48(Proj.=10A,.=2X,„) 




Figure 1. ^pp^ without and with projection for 
Ns = 32 and with projection for Ns = 48. The 
fitting range for rrij-cs is 4 < < < 20. 
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Figure g shows our main result: the comparison 
of the behaviour of the residual mass against Ng 
for different improvement methods. The result 
at small Ng, for example Ng = 16, reveals that 
the projection method does not show any effect. 
In the large-TVs region, the effect of the projec- 
tion of the low-lying eigenvalues on the residual 
mass is, however, remarkable. This is explained 
|5| by simple and qualitative arguments, as fol- 



lows; rrir, 



J dap{a)e 



-aN, 



, where 



a is the square root of eigenvalues of the opera- 
tor A^A and p{a) is the eigenvalue density, which 
grows with a. In the contribution to the residual 
mass of each eigenvalue, there is a competition 
between the density function and the exponen- 
tial suppression. Therefore, at large Ng, the low- 
lying eigenvalues dominate owing to the exponen- 
tial suppression. At small Ng, however, not only 
low-lying but also bulk eigenvalues contribute to 
TOrcs, owing to p{a). 
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Figure 2. Residual mass as a function of Ng with 
and without using the projection method. Filled 
symbols represent our data. We also add data 
from simulations of the CP-PACS and RBC col- 
laborations for comparison. The lines are just to 
guide the eyes. 



5. Summary 

We have studied the residual mass effect of the 
projection method |^ on the restoration of chi- 
ral symmetry for quenched DWF. At large Ng, 
the low-lying eigenvalues of A'^A dominate the 
behaviour of the residual mass. Therefore the de- 
cay rate of the residual mass with Ng can be sim- 
ply controlled by the number of modes projected 
out. The residual mass for smaller Ng, however, 
is controlled by the low-lying and bulk eigenval- 
ues. This is consistent with the arguments in the 
previous section Our numerical study shows 
that the projection method is superior to the im- 
provement of the gauge action at large Ng. 

We have also observed that the quantity 
{J^qP) / {PP) becomes much more stable after 
performing the projection, which may affect also 
other correlation functions to be stable and clean. 

Since the projection method leads to a small 
numerical overhead, we conclude that using Wil- 
son gauge action combined with the projec- 
tion method is competitive with using improved 
gauge actions. To understand the projection 
method further, we are investigating the projec- 
tion method with improved gauge actions [Tl| ]. 
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